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Using the density-matrix renormalization-group, we investigate the critical behavior ol the 
anisotropic Heisenberg chains with spins up to S = 9/2. We show that through the relations 
arising from the conformal invariance and the DMRG technique it is possible to obtain accurate 
finite-size estimates of the conformal anomaly c, the sound velocity v s , the anomalous dimension 
Xbuik, and the surface exponent x s of the anisotropic spin-5 Heisenberg chains with relatively good 
accuracy without fitting parameters. Our results indicate that the entanglement entropy S(L, I a, S) 
of the spin-S 1 Heisenberg chains satisfies the relation S(L, Ia, S) — S(L, Ia, S — 1) = 1/(25 + 1) for 
S > 3/2 in the thermodynamic limit. 
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I. INTRODUCTION 

Conformal invariance plays an important role in the 
study of critical one-dimensional quantum systems. Sev- 
eral constraints appear if we assume that the critical sys- 
tems are conformal invariant ^ In particular, the pos- 
sible classes of critical behavior of the one-dimensional 
quantum systems are indexed by the conformal anomaly 
(or central charge) c as well as the values the anomalous 
dimensions x" ulk of the primary scaling operators O a 

In the mid eighties of last century, it was shown that 
the conformal anomaly c can be extracted from the large- 
ly behavior of the ground state energy E (L). The ground 
state energy of a system of size L behaves as^£ 
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where the constant 5 = 1(4) for the systems under peri- 
odic (open/fixed) boundary conditions, v s is the sound 
velocity, e M is the bulk ground state energy per site, 
and foo is the surface free energy, which vanishes for the 
systems under periodic boundary conditions (PBC). The 
structure of the higher energy states, for a system with 
periodic (open) boundary conditions, are related with 
the anomalous dimensions x^ ulk (surface exponents xf). 
There are a tower of states in the spectrum of the Hamil- 
tonian with energies E^ m , (L) given by^i 
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E« tm ,(L)-E (L) = -^(x + m + m'), (2) 

where m, ml = 0, 1, 2, the constant r\ = 1(2) and x = 
x buik( x< s) f° r the systems with periodic (open) boundary 
conditions. 

The critical behavior of several models were studied 
using the conformal invariance relations above. In par- 
ticular, the spin-1/2 XX Z chain is exactly soluble. For 
this reason, in this case, it is possible to obtain the spec- 
trum of the energies of very large system sizes (solving 



the Bethe ansatz equations). Due to this fact, very ac- 
curate estimates of c, v s , xf ulk and x" were obtained for 
the spin-1/2 XX Z chain4^ However, very few models 
are exactly soluble. And, in general, it is possible obtain 
the eigenspectrum only by exact numerical diagonaliza- 
tions. 

The Lanczos method of exact diagonalization (ED), 
can be used to extract the ground-state energies of the 
Halmiltonians. However, with ED it is possible to con- 
sider small system sizes, since the Hilbert space grows ex- 
ponentially with the system size. Even though, relative 
good estimates of conformal anomaly c and the anoma- 
lous dimension s xZ.n. can be obtained with ED (see, for 
example, Refs. [lfjllll ). 

The density- matrix renormalization-group (DMRG) 12 
is a powerful numerical technique that can be used to 
study very large one-dimensional systems. In principle, 
we can use the DMRG technique and the conformal in- 
variance relations to obtain high accuracy estimates of 
c, x" u[k and x". Since with the DMRG it is possible to 
obtain the ground states energies of large system sizes 
in a controlled way. However, it is not easy to estimate 
the sound velocity v s with the DMRG technique, differ- 
ently of the ED (it is possible to estimate v s with the ED 
by the behavior of the energies in different momentum 
sectors). This is the major reason why the conformal in- 
variance relations above were not explored in the studies 
of critical systems by the DMRG technique.— 

In recent years, this has changed because Calabrese 
and Cardy related the entanglement entropy of one- 
dimensional critical s yste ms with the conformal anomaly 
oi^ (see also Refs. fl5H17l ). As we will see, due to this new 
relation between the confomal anomaly and the entangle- 
ment entropy, it is now possible to obtain the finite-size 
estimates of c, v s , Xbuik , and x s in a systematic way with 
the DMRG technique. The entanglement entropy can be 
denned as following. 18 Consider a one-dimensional sys- 
tem with size L and composed by two subsystems A and 
B of sizes I a and Ib = L — Ia, respectively. The entan- 
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glement entropy is denned as the von Neumann entropy 
S(L,Ia) = —TrpA\ri(pA), associated with the reduced 
density matrix pa- For the critical one-dimensional sys- 
tems the entanglement entropy behaves asi^ 

S(L, l A ) = ^ In sin {^A ) + ci - (1 - V )s b , (3) 

where s;, is the boundary entropy,— c\ is a non-universal 
constant and rj = 1(2) for the systems under periodic 
(open/fixed) boundary conditions. 

Note that we can estimate the conformal anomaly c us- 
ing Eq. [3] without the knowledge of the sound velocity v s . 
Indeed, some authors have plotted S(L, Ia) as function 
of In (— sin f 2 ^ 1 )) in order to estimate c by a numerical 
fit. Here, we use basically the same route that Zhao et 
al. used in Ref. [l9| to extract c in an investigation of the 
XXZ Heisenberg chains with defects. We will estimate c 
by measuring the entanglement entropy of two systems 
with different sizes. With this procedure, we can obtain 
a finite-size estimate of c without a numerical fit. We 
have observed that this procedure gives errors smaller 
than the ones obtained by the fitting procedure (in same 
cases, the error is one order of magnitude smaller). 

It is important to point out that some authors have 
used the DMRG to obtain the critical exponents of sev- 
eral models^ through the asymptotic behavior of static 
correlation functions. This numerical procedure, use 
more CPU time and, usually, provides only rough es- 
timates of the critical exponents. Thus, we see that a 
procedure to obtain accurate estimates of the critical ex- 
ponents (based in the energies behavior) with the DMRG 
is also highly desirable. 

The main aim of this article is to show that if we ex- 
plore all the three conformal invariance relations above 
(Eqs. (l)-(3)) with the DMRG technique, then it is possi- 
ble to obtain accurate finite-size estimates of c, v s , xtuik, 
and x s in a systematic way with the DMRG technique 
without any fitting parameter. Several authors have ex- 
plored some of these relations with the DMRG. However, 
they did not make use of all the three relations (with the 
exception of Ref. l2ll ). We would like to point out here 
that with the multiscale entanglement renormalization 
ansatz (MERA) it is also possible to extract these quan- 
tities with a good accuracy^ 2 - (see also Ref. l23h . 

In order to illustrate the procedure that we use to ex- 
tract c, v s , x^ ulk and x", we consider the anisotropic 
spin- S Heisenberg chains defined as 

H = E ( s Phi + + A*W +1 ) , (4) 

j 

were A = cos(7) is the anisotropy. It is well know that 
this model at the isotropic point A = 1 is gapless (gapful) 
for half-integer (integer) spins^ For this reason, we will 
consider only the half-integer spin cases. Although even 
the models with integer spins become gapless at some 
critical anisotropy^ 



We investigate the model defined above using the 
DMRG 12 method under open boundary conditions 
(OBC) and PBC, keeping up to m=4000 states per block 
in the final sweep. We have done ^6 — 9 sweeps, and the 
discarded weight was typically 1CP 9 — 10~ 12 at that final 
sweep. The dimension of the superblock in the last sweep 
can reach up to 80 millions. In our DMRG procedure the 
center blocks are composed of (2S+1) states. 



II. RESULTS 
A. S=l/2 

We will focus, firstly, in the anisotropic spin-1/2 
Heisenberg chain, which is exactly soluble (for a review 
see Ref. §). For this case, we be able to compare our 
numerical estimates of c, v s , x£ ulk and x" with the ex- 
act results. In order to avoid the logarithmic corrections 
(due to the marginally irrelevant operator) which makes 
the finite-size analysis more complicated, we decided to 
investigate the model at the anisotropic point 7 = 7r/3 
(which corresponds to A = 1/2). Few results for other 
values of 7 were also explored in this work. 

Before presenting the numerical results for the spin-1/2 
chain, let us describe some known results based on ana- 
lytical approaches A 2 ^— The anisotropic spin-1/2 chain 
is critical for — 1 < A < 1 with central charge c = 1, 
the anomalous dimension (surface exponent) associated 
to the lowest eigenenergy in the sector with total spin 
z-component S z = 1 is given by x bu i k = (^s = 2L 7 2 ) 
and the sound velocity is v s = ZLSMll. Very interesting 
to note that for 7 = 7r/3 the anisotropic spin-1/2 Heisen- 
berg chain present a very peculiar ground state. For this 
reason, some analytical expression for the two-point cor- 
relation functions^ and for the reduced density matrix 
were obtained^ Some exact results for the reduced den- 
sity matrix with arbitrary 7 were also obtained by Alba 
et al. in Ref. [3ll 

In Fig. 1(a), we present the entanglement entropy 
S(L,Ia) as function of the length Ia for the anisotropic 
spin-1/2 Heisenberg chain with PBC for a system of size 
L = 100 and 7 = n/3M The dashed line in this figure is 
a fit to our data using Eq. [3] We show S(L, Ia) only for 
l A < L/2 since S(L,L- l A ) = S(L,l A )- The conformal 
anomaly obtained through this fit is c = 1.005, which is 
very close to the exact value c exac t = 1- Estimates of 
the conformal anomaly can also be acquired if we con- 
sider the system with OBC instead of PBC. In Fig. 1(b), 
we show S(L, Ia) as function of Ia for the model with 
OBC, 7 = tt/3 and L = 160. Observe that for the OBC 
case, S{L 1 Ia) exhibits even-odd oscillations as function 
of I a- In this case, for Ia even (odd) the fit to our data 
using Eq. [3] is c = 1.09 (c = 0.77). These even-odd os- 
cillations have been reported in several works *22r— and 
it is expected to decay away from the boundary with a 
power law^ Very interesting to note that even for sys- 
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Figure 1: (Color online). The entanglement entropy S(L,Ia) 
vs. Ia for the anisotropic spin-1/2 Heisenberg with 7 = ir/3. 
The circles are numerical data and the red lines are fits to our 
data using Eq. [3] (see text), (a) S(L, Ia) vs. Ia for a chain of 
size L — 100 and PBC. (b) S(L, Ia) vs. Ia for a chain of size 
L = 160 and OBC. 



tems under PBC, the Renyi entropy also shows even-odd 
oscillations, which have been proposed to have an uni- 
versal scaling l&wM- 

The above route have been used by some authors 
to estimate the conformal anomaly of some criti- 
cal models such as the transverse Ising chain^ the 
isotropic/anisotropic spin-1/2 Heisenberg chahi ) 21 i 36 the 
SU(N) chains*^ and as well as a spin-3/2 fermionic cold 
atoms with attractive interactions^? In order to obtain 
good estimates of c with this procedure we need to con- 
sider large system sizes, mainly for the OBC case, which 
can demand large computational effort in general. Here, 
we use a slightly different route to obtain the finite-size 
estimates of the conformal anomaly c, as mentioned be- 
fore. In particular, for the OBC case, we will see that 
with an extrapolation of the finite-size data, we are able 
to obtain estimates of c and x s with good accuracy. Even 
considering relatively small system sizes, which requires 
a reasonably small computational effort. 

A simply way to extract the conformal anomaly from 
Eq. [31 without any fitting parameter, is considering two 
systems with sizes L and U . Let us assume that these 
two systems are composed of two subsystems of sizes I a = 



Table I: Finite-size estimates of the conformal anomaly c , 
the sound velocity v a , and the anomalous dimension Xbuik 
for the spin-1/2 anisotropic Heisenberg chain with PBC and 
7 = tt/3. We use L' = L + 16 in Eq. [5] 



L 


C PBC 


v a 


Xbuik 


16 


0.99596 


1.30675 


0.33314 


32 


0.99879 


1.30103 


0.33324 


48 


0.99942 


1.29994 


0.33328 


64 


0.99966 


1.29955 


0.33330 


80 


0.99977 


1.29937 


0.33331 


96 




1.29927 


0.33332 


exact 


1 


1.29903 


0.33333 



L/2 and l A = L' /2, respectively. Thus, from Eq. [31 we 
see that we can estimate c by 



c{L,L') = 3r] 



S(L,L/2)- S{L',L'/2) 3//AS 
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(5) 



Note that previously work a 19 ' 38 also use the incre- 
ment of the entropy, AS, to extract the conformal 
anomaly of the spin-1/2 XXZ Heisenberg model with de- 
fects/impurities. Lauchli and Kollath 3 ^ also used the in- 
crement of entropy in order to locate the quantum critical 
point of the Bose-Hubbard chain. 

In Table I, we show the finite-size estimates of the con- 
formal anomaly, obtained by Eq. [5J for the anisotropic 
spin-1/2 Heisenberg chain with PBC and 7 = 7r/3. Sim- 
ilar results were also found for 7 = ir/6 and 7 = ir/8. 
For comparison purpose, we also present in this table the 
exact values of the conformal anomaly c, the sound ve- 
locity v s , and the anomalous dimension Xbuik- Note that 
by using relatively small system sizes, we are able to esti- 
mate c with a small error (~2xl(T 4 ) without any fitting 
parameter. 

Once the conformal anomaly c is obtained, we can use 
Eq. [1] to extract the sound velocity v s . In particular, 
for the PBC case the finite-size estimate of the sound 
velocity is obtained by^ 



v s (L) = — (Le c 
ire 



E (L)). 



(6) 



Note that in this equation, we also need the bulk ground- 
state energy per site eoo. However, this is not a problem 
since with the DMRG technique it is possible to obtain 
eoo with a high accuracy. 

The finite-size estimates of v s obtained using Eq. [6] 
are also presented in Table I. As observed in this table, 
accurate results are also acquired for v s . As we already 
mentioned in the introduction, until very recently, it was 
not possible to extract the sound velocity v s , based only 
in the large-L behavior of the ground state energy with 
the DMRG technique. But, as we have observed, with 
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Table II: Extrapolated and finite-size estimates of the con- 
formal anomaly c OBC and the surface exponent x a for the 
anisotropic spin-1/2 Heisenberg chain with OBC and 7 = 7r/3. 
We use L' = L + 20 in Eq. [S] The extrapolated values were 
obtained by a numerical fit (see text). 



L 


c OBC 




60 


1.0859 


0.64298 


80 


1.0727 0.64809 


100 


1.0633 


0.65131 


120 


1.0568 


0.65353 


140 


1.0518 


0.65517 


160 


1.0470 


0.65642 


00 


1.004 


0.66658 


exact 


1 


0.66666 



As we see in this table, the estimate of c obtained 
with this procedure gives an absolute error of about 10~ 3 , 
which is quite good. 

In a similar way, we also use Eq. [7]to obtain the surface 
exponent x s . In Table II, we also present the finite-size 
estimates of x s for the anisotropic spin-1/2 Heisenberg 
chain with OBC and 7 = 7r/3. The extrapolated value of 
surface exponent x s was obtained by a numerical fit, as 
we did for the conformal anomaly. We assume that the 
surface exponent behaves as 

x s (L) =x s +a/L + b/L 2 . (9) 

As observed in the table, the extrapolated value of x s is 
also very close to the exact one. 

B. S>l/2 



the use of the Eqs. Q] and we can estimate v s , with 
good accuracy, considering relatively small system sizes. 

Finally, using Eq. we extract the finite-size estimates 
of the anomalous dimension/surface exponent by 



x(L) 



rjL (Ei(L) - E (L)) 
2nv s 



(7) 



where E\ (L) is the ground state energy in the sector with 
total spin z-component S z = 1. In Table 1, we also show 
the finite-size estimates of Xbulk for the anisotropic spin- 
1/2 Heisenberg chain under PBC and 7 = 7r/3. Note 
again that accurate results are obtained. 

From the above results, we learn that it is possible 
to obtain good estimates of c, v s and x^uik, without 
any fitting parameters (or extrapolations), considering 
relatively small system sizes under PBC. On the other 
hand, for the OBC case is necessary to extrapolate the 
finite-size estimates in order to obtain accurate results, 
as shown below. 

Let us now estimate the conformal anomaly consider- 
ing the model with OBC (Eq. [5] with 77 = 2). In Table 
II, we present the finite-size estimates of the conformal 
anomaly c for the anisotropic spin-1/2 Heisenberg chain 
with OBC and 7 = 7r/3. As we observed in this table, 
the finite-size estimates of c are less accurate than 
the ones found for the model with PBC, even consider- 
ing larger systems for the OBC case. We can improve 
the estimate of c OBC extrapolating the finite-size data 
to the infinite lattice value. We assume that conformal 
anomaly behaves as 



„OBC 



(L) =c + a/L + b/L 2 



(8) 



The extrapolated values in Table II, were obtained from 
a fit to our data with this equation (we use only the last 
four points of our data in the numerical fit). Note that 
the extrapolated estimate is quite close to the exact value. 



The results obtained above for the spin-1/2 chain are 
all known, and we studied the spin-1/2 case for the pur- 
pose of comparison/benchmark. However, that study 
was highly important, because it established that the pro- 
cedure used by us provides accurate results. Now, let us 
consider the cases where S > 1/2. 

Before presenting the estimates of c, Xbuik, and x s . We 
will present an interesting behavior of the entanglement 
entropy S(L, I a, S) for the spin-S 1 Heisenberg chains with 
spin S > 1/2. We observed what appears to be an uni- 
versal behavior of S(L, I a, S) for S > 1/2. Our numerical 
data indicate that the entanglement entropy of two sys- 
tems with spins S and S — 1 respectively, with S > 3/2 
are related by following the relation 



S(L,l A ,S)-S(L,l A ,S-l) 



(25-1) 



ai, 



(10) 



where a\ is a very small term. Note that this is equivalent 
to say that the non-universal constant C\{S) (see Eq. [3]) 
satisfies the relation C\{S) — Ci(S — 1) = 



a 1 . 



(25-1) 

We have observed that constant at decrease with the 
size of the systems. We were not able to proof that a\ is 
zero in the thermodynamic limit, although our numerical 
results suggest that. In particular, for a system with PBC 
(OBC), 7 = tt/3, and size L — 80 (L=160) we found 
that a± ~ 10 _3 (ai ~ 1CP 2 ). In the discussion below we 
neglect this term. Due to the Eq. [TUl we can express the 
entanglement entropy S(L, I a, S) of a systems with spin 
S > 3/2, in terms of S(L, Ia, 3/2) by the equation 



(S-5/2) 



S(L,l A ,S)=S(L,l A ,3/2)+ J2 



1 



H 2 ( 5 -i)- 1- 
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In order to check the valid of the Eq. [Tl] we investigate 
the following function 

(S-5/2) 

W(L,I a ,S) = Sdmrg(L,I a ,S)- 



1 
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Figure 2: (Color online). The rescaled entanglement entropy 
W(L, Ia,S) vs. I a for the anisotropic spin-S Heisenberg for 
some values of S (see legend), (a) W(L,Ia, S) vs. Ia for 
chains of size L = 80, 7 = tt/3 and PBC. (b) W(L,Za, S) vs. 
Za for chains of size L = 32, PBC and some values of 7 (see 
legend). The simbols refer to the values of the spins (we used 
the same simbols of figure (a)), (c) W(L,Ia, S) vs. Ia for 
chains of size L = 160 and OBC. 



and we set W(L,l A ,3/2) = S DMRG {LJ A ,3/2). 

Let us focus, firstly, in the PBC case. In Fig. |5[a), we 
present the rescaled entanglement entropy W(L,l A ,S) 
as function of Ia for some values of spins for systems 
with PBC, L = 80 and 7 = tt/3. As observed in this 
figure, all curves collapse onto a single universal scaling 
curve. Actually, there is a small difference around 10~ 3 
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Figure 3: (Color online). Some results for the anisotropic 
spin-S Heisenberg with OBC, 7 = 7r/3 and some values of S 
(see legend), (a) The entanglement entropy S(L,Ia) vs. Ia 
for chains of size L = 160. (b) The dimer parameter D(j) vs. 
j for chains of size L = 80. We present D(j) only for distances 
j<40. 



between these curves (see Figl^b)). Similar results are 
also observed for other values of 7 for systems of size L — 
32, as shown in Fig. [5{b). Note that as 7 decrease the 
differences between the curves increase. This is expected 
since finite-size effects become stronger as 7 — > (or A — > 
1 ) . These results strongly indicate that Eq. [11] is valid 
for any values of spin S > 3/2 and 7. Our numerical 
data (not shown) support also that Sdmrg(L, l A ,3/2) = 
Sdmrg(L, Ia, l/2) + ll/18 for the PBC case. We will see 
below that if we consider the systems with OBC, the Eq. 
[□still holds true for S > 3/2. 

Now, let us consider the OBC case. In Fig. [3ja) we 
show S(L, Ia) for the anisotropic spin-S Heisenberg with 
OBC, L = 160, 7 = 7r/3 and some values of spins S. 
Note that differently of the spin-1/2 case with OBC, for 
spins S > 1/2 the entanglement entropy S(L,Ia) does 
not present the strong odd-even oscillations. Laflorcn- 
cie et al. in Ref. |33j argue that the oscillating term 
in S(L, Ia) is universal, and is connected with the an- 
tiferromagnetic nature of the Hamiltonian (however, at 
the present moment, there is no analytical derivation of 
this statement). In particular, the alternating term in 
S(L, l A ) appears to be proportional to the alternating 
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term in the energy density 33 

h(j) =< ( S | S j +1 + a y j+1 + a S | S | +1 ) > . (is) 

In order to understand the absence of the strong odd- 
even oscillations in entanglement entropy S(L, l A ) for 
S > 1/2, we study the behavior of the "dimer parameter" 
D(j) = h(j + 1) — h(j). It is expected that D(j) decay 
with a power law^ 3 - In Fig. [3Jb), we present the dimer 
parameter D(j) for chains of size L = 804i We roughly 
estimate the errors of D(j) around ~ 10~ 6 , since we find 
numerically that (Eq — Ylj h(J)) ~ 10~ 6 . As observed in 
this figure, the magnitude of the oscillations of h(i) for 
S > 1/2 are very small. These results, indeed, corrobo- 
rate that the alternating terms in S(L, l A ) and h(J,A) are 
connected. However, due to the fact that the values of 
D(i) is so small for S > 1/2, we were not able to check 
that the alternate parts of S(L, I a) and h(l A ) decay with 
the same exponent, as expected by Laflorencie et al££ 

As in the PBC case, we also observed that the curves of 
the rescaled entanglement entropies W(L, I a, S) of chains 
with S > 3/2 and OBC collapse onto a single universal 
scaling curve, as we can see in Fig. [He). 

Finally, let us present our estimates of c PB , Xbuik and 
x s for several values of spins S. In Table III, a summary 
of our results is provided (acquired following the pro- 
cedure explained in the subsection II. A). The finite-size 
estimates of c PB and Xbuik were obtained considering 
systems of sizes L = 80 or L = 96 with PBC, while the 
extrapolated values of c OBC and x s were obtained con- 
sidering the model with OBC and sizes up to L = 180. 
We also show in this table the ratios of the dimensions 
x s / Xbuik- Based in our benchmark results of the spin-1/2 
case, we believe the errors of the quantities presented in 
table III are smaller than 10~ 3 . 

Note that it is expected that the critical behavior of 
the Heisenberg chains with half-integer spins belongs 
to the same class of the Gaussian model s 10 ' 42 In par- 
ticular, numerica l 8 ' 10 ' 43 and analytica l 25 ' 28 ' 42 techniques 
show that c = 1 for the Heisenberg chains with half- 
integer spins. Moreover, it is expected that anomalous 
dimension of the anisotropic Heisenberg chains depend 
of the anisotropy and of the value of the spin 5*i^ i 10 ' 28 ' 44 

As observed in Table III, our estimates of c are in 
agreement with the expected value of c = 1. More- 
over, our results for 5 = 3/2 are in perfect agreement 
with the ones found by Alcaraz and Moreoii - In partic- 
ular, for 7 = 7r/3 they found the following extrapolated 
values: c = 1.08, Xbuik — 0.098 and x s = 0.198. Our re- 
sults are also consist with the Alcaraz-Moreo's conjecture 
x s = Ixbuik- This conjecture was proposed based in the 
exact diagonalization calculations of small system sizes. 
Besides that, Alcaraz and Moreoi - reported results only 
for spins up to S = 2 for the antifcrromagnetic region 



< A < 14 5 - Here, considering larger system sizes (giv- 
ing better estimates) and larger values of spins we show 
that this conjecture holds. 

Table III: Finite-size estimates of the conformal anomaly 
c PBC and the anomalous dimension Xbuik for the anisotropic 
spin-S Heisenberg chain with PBC and 7 = 7r/3 obtained with 
L — 80-96. The extrapolated values of the conformal anomaly 
c OBC and the surface exponent x s are also presented for the 
same model/coupling with OBC. 



O C C Xbuik X s 

1/2 0.9997 1.004 0.3333 0.6665 1.999 

3/2 0.9995 1.002 0.09918 0.1984 2.000 

5/2 0.9993 0.999 0.0572 0.1143 1.998 

7/2 0.9989 1.002 0.0403 0.0806 2.000 

9/2 0.9995 1.001 0.0311 0.0623 2.003 



III. CONCLUSION 

In this article, we present a simply procedure to obtain 
accurate estimates of the conformal anomaly c, sound ve- 
locity v s , anomalous dimension Xbuik and the surface ex- 
ponent x s using the conformal invariance relations and 
the density-matrix renormalization-group technique. In 
order to illustrate the procedure we use to get these quan- 
tities, we investigate the anisotropic spin-S Heisenberg 
chains with periodic/open boundary conditions. 

Our results for the model with spin S=l/2 were com- 
pared with the exact results. For the spin-1/2 case, we 
found that the procedure that we used to estimate c, v s , 
Xbuik and x s gives accurate estimates with errors smaller 
than 10 -3 . We also present accurate results for the model 
with spins S = 3/2, 5/2, 7/2 and 9/2, and we confirm 
the Alcaraz-Moreo's conjecture a; s = 2xbuik- 

Our numerical results also support that the entangle- 
ment entropy of the spin-S 1 Heisenberg chains satisfies 
the relation S(L,l A ,S) - S(L,l A ,S- 1) = 1/(25*+ 1) 
for S > 3/2 in the thermodynamic limit. We also veri- 
fied that the alternate terms of S(L, Ia) and h(l A ) seem 
to actually be connected, as suggested by Laflorencie et 
alJ£ However, we were unable to verify that both terms 
decay with the same universal exponent. In this vein, 
it is interesting to observe the universal oscillatory be- 
havior of the Renyi entropy predicted by Calabrese and 
co-authors in Ref. This study is in progress and the 
results will be present elsewhere. 
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